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Abstract 

We present new periodic, kink-like and soliton-like travelling wave so- 
lutions to the hyperbolic generalization of Burgers equation. To obtain 
them, we employ the classical and generalized symmetry methods and the 
ansatz-based approach 

1 Introduction 

Actually there are no analytical methods enabling to solve an arbitrary initial or 
boundary value problems for nonlinear PDEs that are not completely integrable. 
The majority of evolutionary PDEs used to simulate non-linear transport phe- 
nomena is not completely integrable. Yet it is of common knowledge that under 
certain conditions coherent structures formation take place in open dissipative 
systems, being simulated by PDEs, that are not completely integrable. The 
analytical description of the coherent structures within the non-integrable mod- 
els is of great interest, since analytical form of solution is more preferable in 
analyzing structures formation and evolution and, besides, they are widely used 
as a starting point for various asymptotic methods, for testing the numerical 
schemes and facilitating the stability analysis. 

To obtain exact solutions to non-linear PDEs, various methods are used. 
Besides the classical symmetry reduction scheme (see e.g. which is not 

effective in obtaining solutions with the given properties, there are employed 
techniques based on choosing a proper transformation (or ansatz), turning the 
problem of finding out exact solutions to the algebraic one EJ E] • In our 
previous study [5] we employed a certain dialect of the the direct algebraic 
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method to the hyperbolic generalization of Burgers equation (GBE). Yet in 
the following publications [3] we showed that within this methodology it 
is impossible to obtain the the solitary wave solution of GBE, occurring to 
exist for certain values of the parameters [7J- This served us as a motivation 
for the developing the effective methods of approximations [7| |S] of solutions 
describing the coherent structures. Actually we return to the subject of the 
analytical description of the coherent structures within the GBE model, using for 
this purpose quite different technique, namely, a generalized symmetry and the 
combination of the classical symmetry reduction with an ansatz-based method 

M 

The structure of the study is following. In section 2 we show how it is pos- 
sible, combining the classical reduction with the Hirota-like ansatz, to linearize 
equation describing a set of travelling wave solution to GBE and obtain on this 
basis new periodic, kink-like and soliton-like solution. In the following section 
we show how it is possible to gain the linearization of the reduced GBE, using 
the Painleve test. In section 4, using the idea of conditional symmetry, we obtain 
a wide class of exact solutions to GBE, containing the functional arbitrariness. 
In section 5 we give some examples of such solutions. In the final part of the 
study we make some conclusions and general remarks. 



2 Solutions to GBE obtained via the classical 
symmetry reduction and Hirota-like ansatz 

Let us consider a hyperbolic generalization of Burgers equation |1U1 |S] : 

Tu tt + Auu x + Bu t + Hu x - ku xx = f(u) — A (u - mi)(u - m 2 )(u - m 3 ). (1) 

We are going to look for the wave patterns among the set of travelling wave 
(TW) solutions, i.e. solutions invariant with respect to translations generators 
-^j and t^:. So it is quite natural to employ the classical reduction scheme 
12, and consider an over-determined system containing action of the generator 
X = Sr — u ■§- on the solutions of GBE: 

at ^ ox 

Tu tt + Auu x + Bu t + Hu x - ku xx = f(u), (2) 

X [u(t, x) — u] = ut — [iu x = 0. (3) 

The classical theorems assure that system ©-Q is compatible and has 
non-trivial solutions. In fact, solving equation J3J we obtain an ansatz 

u(t,a:) = Z7(0, £ = x + »t. (4) 
Inserting ansatz (@J into J2J, we obtain the following ODE: 

(tu 2 - K ) U" (0+(H+B fx) U'(t)+AU(£) U'(0 = A (t/(£)- mi )(tf(£)- m2 )([/(£)- m3 ). 

(5) 
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The problem is that, generally speaking, equation JSJ is non-integrable, so, 
employing the classical reduction scheme we are not able to obtain analytical 
description to wave patterns. 

In the following we use this scheme together with the Hirota-like ansatz 

*'(£) 

U® = , (6) 

and show that under certain conditions such combination leads to a linear ODE. 
Inserting © into (JjJJ, we obtain the equation 

*(£) 2 [A mim 2 m 3 *(£) - A (m 2 m 3 + mim 2 + mim 3 )*'(f) + (H + 
-A (m x + m 2 + ma) + B/*)] + [*'(0] 3 {A-\-2k + 2/j 2 t) = 0. 

(7) 

One easily gets convinced by the direct inspection that the following statement 
holds: 

Lemma 1 Ansatz @) leads to the linear ODE 

A [mi 77127713^ — (777,27713 + 7711777,2 + mim3)^f']+(H+B fj,)^/" +(fi 2 r— n)^'" = 0, 

(8) 

provided that following conditions are fulfilled: 

A + 3k - 3/j 2 t = 0, (9) 

H - A (mi + 777,2 + 77i 3 ) + Bfj, = 0, (10) 
A - A - 2k + 2/j 2 t = 0. (11) 

Equation (jHJ is a third-order ordinary linear differential equation. Its solu- 
tions depend on the roots {cfc} fc=1 2 3 of corresponding characteristic equation. 
In case when the conditions ijjfjl— Qlljl- arc fulfilled, the characteristic equation is 
as follows: 

0~ — (mi + 777,2 + 77l3) (7 2 + (777,2 77i3 + 777,1 77i2 + 777l 7773) °~ ~ m l m 2 7773 = 0. (12) 

The roots of equation l|12|) coincide with the numbers mi, 7772 and 777,3, being 
the roots of equation f(u) — 0. Below we consider four distinct cases. 

Case I. For mi 7^ m 2 7^ 7773 7^ mi the general solution of © takes on the 
form 

= e mi « Cl + e m2 «c 2 + e m3 «c 3 . 
In the following we put c% = 1. With this assumption we obtain the expression 

777ie mi ? + c 2 m 2 e m2 ^ + c 3 m 3 e m3 ^ , . 

u(t, x) — 1 - t t , F = x + fit. (13) 
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Figure 1: An example of TW solution described by the formula l|13fl 



For positive c 2 and c 3 formula (| 1 31) describes a kink-like solution (see fig ^| . 
Case II: For wi\ 7^ m 2 = m 3 

*(0 = e mi « + e" l2 «(c 2 + c 3 

and using the formula 10 we get the solution 

mie mi5 + C3e m2 «(c2m2 + C3 + TO2C 3 £) 

= 1 T, Tx ■ (14) 

e mi ? + e m2 «(c 2 + c 3 £) 

Asymptotic analysis shows, that solution (|14|l tends to either m\ or to m 2 
as £ — > ± 00. The condition of non-singularity of the solution reads as follows: 
V(£i?: e miC + e m2? (c 2 + c 3 £) > or, in other words, 

e( m ™« >-(c 2 + c 3 0- (15) 

It is possible when mi — m 2 is positive, and c 3 is negative or when m 2 — mi, and 
c 3 are positive. In fig. |2 the first possibility is illustrated. Note that the second 
one is symmetric. It is evident from the analysis of fig. [21 that condition (|15|l is 
satisfied in case when — c 2 < c cr . One can easily verify that c cr is expressed by 
the formula 



-A[c 3 



■ e 



m 2 l 



providing that A = — - — In — — — . 

1 mi— ni2 m2— mi 

Case III: For mi = ?ti 2 = m 3 = m 
and using (j^J we get the solution 

u{t ' x) = m+ c 3 + a^Uy (16) 

Formula describes a a solitary wave with " heavy" tails (fig. 0J| , providing 
that inequality c| — 4 c 3 < holds. 
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Figure 2: Case II: nonsingular solutions. Solid lines correspond to the plots of 
functions /i(£) = e^™ 1- ™ 2 ^ and /2(C) = ^( c 2 + ca£). Dashed line is a tangent 
to being parallel to / 2 (£) 




Figure 3: An example of TW solution described by the formula 114fl 
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Figure 4: An example of TW solution described by the formula i|16fl 
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Case IV: For mi, m 2 = a ± /3z, and real 7713, we get the solution 

,~ = c 3 m 3 e m ^ + 2e^ [acos(/?Q - /?sin(/3fl] 

c 3 e™3? + 2e««cos(/3£) 1 J 

Let us notice, that the only possibility for the solution p7|) to be non-singular 
is simultaneous fulfillment of the conditions m 3 = a and | C3 1 > 2. In this case 
we get a periodic solution (fig. |5J| . 

3 Painleve analysis of GBE 

Now we are going to show how the conditions on constants similar to l|9")l- Hll[) 
can be obtained using the Painleve test 1 1 1 j . Since we concentrate here merely 
upon the set of TW solutions, let's start from the equation J3l- 

In the first step we take [/(£) = U(x + fit) = a[0]£ p and put it into the 
leading terms of JSJ). This way we obtain the relation: 

{p - p 2 )(p 2 r - k) + Ap£ p+1 a[0] + \£ 2p+2 a[0} 2 = 0. (18) 

In order that the equation be balanced, we must put 2(p + 1) = p+1, so p = —1. 

Next we take £/(£) = X)to a M? l_1 - Equating to zero coefficients of each 
power of £, we get the recurrence for a[i\. The first two expressions are as 

follows: 

k\ a[0] = A ± ^/A 2 - 8A (p 2 T - n) (19) 

(A - 3Aa[0])a[l] = -(H + B\i + A (mi + m 2 + m 3 )a[0]) (20) 

Since in the Painleve test being applied to a second order equation one of a[i] 
should be " free" , then we can put 



A - 3Aa[0] = 

H + Bfi + A (mi + m 2 + m 3 )a[0] = 0. 



(21) 



Solving (|21|l wc obtain: 



A = — , A — —6s, 



/i 2 t — n ' mi + ni2 + m 3 

If we additionally put a[0] = 1 then we get precisely the conditions for which 
equation Q linearizes. 

4 Exact solutions of GBE associated with the 
conditional symmetry 

Now let us consider an over-determined system 

tu u - ku xx +uu x +u t = f(u), (22) 
X [u] = vu t + u x = $(w), (23) 
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Figure 5: Periodic TW solution to J5J) described by the formula ijTTj) 



where 

X = vd t +d x + $(w)<9„ (24) 

is a generalized (conditional) symmetry generator containing an unknown func- 
tion <E>(u). 

In accordance with |12M13llT4*| . system l|22H23[) will be compatible if it admits 
the second prolongation of the operator X: 

d d d d d 
X 2 = X + ( xx +( x +Q t + ( tx (25) 

OU X x OU t OU x OUtt OUtx 

where 

Ct = A $("), ( x =D x $(u), (tt = Dt(t, (tx=D t ( x , ( XX = D X ( X 

and D t , D x are total derivatives with respect to corresponding variables. Ap- 
plying the operator (|2*5ll to equation (|2*2"|> . we get: 

Ct+r Ctt + Hu) u x + u ( x - k Cxx - f(u) $(u)k=o = 0, (26) 

where Ei = means that equation is considered on the manifold, i.e. when 
equations (|22|) ~ (|23[) and all their differential consequences are taken into ac- 
count. 

A passage to the manifold defined by equations (|2*2*|l - lf2l^l requires some 
comments. Equation (|23|) and its differential consequences give us the following 
conditions: 

u x = §{u)-vu t , (27) 
Uxx = $(u) $(u) - 2 v <S>(u) u t + v 2 u tt . (28) 

Using these formulae and equation (I22|) . we can also exclude u u - In fact, insert- 
ing lt77|) - (|2*5|) into we obtain the equation 

(r - kv 2 ) u tt = f{u) + (uu - 1 - 2/£i/$(t()j u t + $(u) (k$(u) - uj . (29) 

Here we have two possibilities, depending on whether (t — k^ 2 ) is zero or not. 
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If t — nv 2 7^ 0, then, returning to Q26[l and taking account of l|27[) - i|28[) and 
(|29() . we obtain a first order PDE. Applying the standard procedure of splitting 
|12l 1131 HI] and solving the system of determining equations we conclude that 
$[it] = and therefore X is a classical symmetry generator. 

Now let us assume that r — kv 2 — 0. Using the formulae (|T7Jl - lj25|) . we 
can exclude from lf^|l merely the spatial derivatives u x and u xx . The splitting 
procedure gives us in this case the following system: 



(30) 
(31) 



k j$(u)$(u) + $(u) | = $(u) + u$(u) - f(u) 
2kv j$(u)<3>(w) + $(u) | = (vu- l)$(u) + v$(u) 

Integrating once equations l|37)|l - (|3"T)) . we get the system 

k$(w)4>(w) = u$(u) - f(u) - Cl, 
2ki/$(m)$(«) = (vu- 1) $(tt) + C 2 . 



(32) 
(33) 



Multiplying i|32[l by 2v and extracting the resulting equation from (|33|l , we get 
the relation between the functions <!>(«,) and f(u): 



$(w) 



2v [f(u)+Ci]+C 2 

V u + 1 



(34) 



Let us note that function $ (u) is defined by the Abel-type equation (|33|) . 
Passing to the new independent variable ( — (vu — 1) C2I (2kv 2 ) , we can 
bring it to the canonical form |15| : 



*(C) 



aC*(C) + l, 



2 k i/ 2 



A passage to the new variables 



leads to the Riccati equation 



e = $--* 2 (0 



d£ 2 ^ 



Solution to this equation is expressed by the following formula |15| : 



* (0 = ve 



AJ_ 



+ BYi 



(35) 



where A, B are constant parameters while J m (z), Y m (z) are the Bessel func- 
tions. 

So, in accordance with the formulae l|35|l and l|34l) . f(u) is expressed in a 
complicated way by special functions. In order to obtain an explicit description 
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to f(u), we set C\ = C2 = 0. In this case system (l32H - (|33fl has the following 
solution: 

f { u) = hQ+ (uh - I ^)u-^ + l-u^ (36) 

where ho is an arbitrary constant. Now we are going to formulate the main 
result of this section. 

Theorem 1 Suppose that f(u) and are given by the formulae ^3b]) and 

{Sty respectively. Then 

• system l^22\) - ^2c!\) is compatible: 

• every solution of $23\) satisfies 



5 Examples of conditionally-invariant solutions 
of GBE 

For v = r = k = 1 equation (|23[1 takes the form 

u t + u x = i (u 2 - 2u + 8h Q ) . (38) 

General solution of this equation depends on whether A = 1 — 8 ho is positive 
or not. For ho < 1/8 solution of Q38JI is as follows: 

x ^ = J t ^ r ' (39) 

G(u)) e 2 - 1 

where 

ui = l + VA, u 2 = 1-VA, 

G(-) is an arbitrary function of u = x — t. Putting ho = —1 and G(ui) — —e r ^ 
we obtain the formula 

V ' ' 1 + exp[3i/2 + r(w)] v ' 

For ho = — 1 and G(w) = e r ^ we get the solution 

v ' y l-eap[3t/2 + r(w)] V ' 



If /to > 1/8 then solution to l|38l) is as follows: 

I fit 

u(t, x) = 1 + j3 arctg 



G(io) 



(42) 
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Figure 6: Temporal evolution of solution described by formula 1401) in case when 
r(w) = sin[2.25w] 
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Figure 7: Temporal evolution of solution described by formula H40|) in case when 

r(w) = -c; 2 



where j3 = •v/|l — 8 ho\. This solution is always singular. 
If ho = 1/8 then solution to Ip^jl is as follows: 

u(t, x) = 1 + 1 t , (43) 

This solution is also singular. 

Let us give examples of solutions corresponding to formulae (|40|1 and (|41|) . 
Thus, inserting r(w) = sin [2.25 w] into equation H41[>. we obtain an oscillating 
kink-like solution, shown in fig. [fj] For T(lu) = —lj 2 this solution produces a 
"dark" soliton with a growing support (fig'EI)- 

In contrast to (|40|l . solution (|41|l is always singular. For r(uj) = — 3.75 lu 2 + 5 
its evolution is shown in fig. |SJ in which we see how an initial localized wave 
pack grows in amplitude and in a finite time gives rise to a blow-up regime. 



6 Concluding remarks 

In this study new solutions to GBE, describing wave patterns formation and 
evolution have been obtained by means of classical and generalized symmetry 
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methods and ansatz-based approach. Most of these solutions cannot be ob- 
tained within the dialects of direct algebraic methods, previously applied to 
this equation [Hj. 

Let us discuss the solutions obtained in section 2. In the case when mi 7^ 
ni2 7^ 7713 7^ mi we deal with solutions described by a rational combination of 
the exponential functions which, generally speaking, do not coincide with the 
powers of a single exponential function as it was the case in . In the following 
cases considered in section 2, solutions are described as rational combinations 
of exponential, polynomial and trigonometric functions and such combinations 
were never applied to GBE before. 

Solutions obtained in section 4 contain arbitrary functions, depending on the 
linear combination of spatio-temporal variables. From one hand, it allows, by 
a proper choice of these functions, to construct a variety of different solutions, 
such as those presented in section 5. From the other hand, appearance of an 
arbitrary function can be employed to describe a sufficiently wide family of 
an initial and (or) boundary value problems. Obtaining of such solutions has 
become possible due to employment of conditional symmetry. 

The results presented in section 2 suggest that the possibilities of obtaining 
new solutions for GBE within the ansatz-based methods are not exhausted. It 
should be noted, however, that effectiveness of these methods grows in the situa- 
tion when they are combined with another methods such as symmetry reduction 
or the Painleve test. As it was shown in section 4, application of the methods 
based on conditional symmetry is also very promising. Let us note in conclu- 
sion, that for GBE the most simple conditional symmetry has been obtained 
as yet. Our preliminary investigations show, that there are another conditional 
symmetry operators admitted by GBE, which seem to be very promising from 
the point of view of obtaining new sets of solutions. 
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